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I ' Abstract 
5-1 ■ 

An 'effective' quasi-local energy expression, motivated by the (relativistically 
corrected) Newtonian theory, is introduced in exact GR as the volume integral of 
^ I all the source terms in the field equation for the Newtonian potential in static space- 

rs ■ times. In particular, we exhibit a new post-Newtonian correction in the source term 

in the field equation for the Newtonian gravitational potential. In asymptotically 
flat spacetimes this expression tends to the ADM energy at spatial infinity as a 
monotonically decreasing set function. We prove its positivity in spherically sym- 
^3 ' metric spacetimes under certain energy conditions, and that its vanishing character- 

izes flatness. We argue that any physically acceptable quasi-local energy expression 
should behave qualitatively like this 'effective' energy expression in this limit. 
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b 1 Introduction 



In non-gravitational classical field theories on fiat Minkowski space the energy-momentum 
distribution of the matter fields is described by the symmetric energy-momentum tensor 
Tab satisfying the dominant energy condition. Then the quasi-local energy of the matter 
fields ElD^K""] with respect to some constant future pointing timelike unit vector field 
K'^ (i.e. a time translational Killing vector of the fiat spacetime) is defined to be the 
integral jj^ KaT°'''tbdT, on the compact domain D G with boundary S := dD in some 
spacelike hypersurface S. Here is the future pointing unit normal to E and dE is the 
natural volume element. As a consequence of the dominant energy condition this is not 
only positive definite, but also is a monotonically increasing set function: if Di C -D2 
then E[Di,K"'] < E[D2,K"']. This implies, in particular, that in asymptotically fiat 
configurations (when E extends to spatial infinity and the total energy -E[E, K"-] is finite) 
E[D, K""] tends to the total energy from below as D is blown up to exhaust E. 
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Although in general relativity there is no well defined energy-momentum density of the 
gravitational 'field', in asymptotically flat conflgurations its total (ADM) energy could 
be deflned, and one of the greatest successes of classical general relativity in the last 
third of the 20th century is certainly the proof by Schoen and Yau [IJ that the total 
gravitational energy is strictly positive deflnite. The logic of one of its simplest proofs, 
due to Witten [2] (and simplifled and corrected by Nester [3j), is that we can rewrite the 
total energy as an integral of some expression (the so-called Sparling form [4J, see also 
[5j) on a spacelike hypersurface, and by Witten's gauge condition the integrand could be 
ensured to be pointwise strictly positive deflnite. Thus the negative deflnite part of the 
Sparling form in the integrand is a pure gauge term. 

These results may yield the view that if the gravitational mass could be deflned at 
the quasi-local level, then it would have to be not only positive deflnite, but also that in 
asymptotically flat spacetimes it would have to tend to the ADM mass as an increasing 
set function. (Note that while we can compare the quasi-local masses, i.e. scalars, we 
can compare the quasi-local energies, i.e. components of four-vectors on different 2- 
surfaces only in the presence of some extra structure, e.g. for spherically symmetric or 
stationary systems.) In fact, several speciflc quasi-local mass expressions exist which 
satisfy these requirements (viz. the Bartnik [6] mass, the mass built from the Dougan- 
Mason [7J energy-momenta, and the Misner-Sharp energy [8j for spherically symmetric 
conflgurations), or at least the second of these in certain special spacetime conflgurations 
(the Hawking [H] or the Geroch energies |1U] and the Penrose mass pjj). However, 
there are other constructions (e.g. the Brown- York expressions [12], the Epp [13], the 
Kijowski-Liu-Yau [T^ [15] and the Wang- Yau energies [16j ) which tend to the ADM 
energy as decreasing set functions. (For a review and a more detailed discussion of the 
various quasi- local energy constructions and the extended literature, see e.g. [IT].) These 
different monotonicity properties of the quasi-local mass/energy expressions generated 
some debate in the relativity community: whether a physically reasonable quasi-local 
energy expression should be monotonically increasing or decreasing near spatial inflnity as 
it tends to the ADM energy. In fact, near spatial inflnity the matter and the radiation can 
be neglected compared to the negative definite Newtonian gravitational binding energy. 
Therefore, increasing the domain of integration we have more and more negative definite 
contribution to the total energy, which therefore must be a decreasing set function (see 
also [18]). 

Since near spatial inflnity the dynamics of the flelds and the gravitation dies off 
rapidly (typically as 1/r^), in the analysis of the asymptotic behaviour of any quasi- 
local energy/mass expression it is natural to consider the spacetime to be static in the 
flrst approximation. The advantage of the existence of a static Killing vector is that it 
provides a geometrically preferred notion of time and a preferred foliation of spacetime 
by a geometrically distinguished family of extrinsically flat spacelike hypersurfaces. In 
fact, the Newtonian limit of general relativity is deflned in this way in [19j. Identifying 
the Newtonian potential with the logarithm of the length of the Killing fleld, we show 
that the exact fleld equations for this take the form of a Poisson equation in which 
the source term contains not only the familiar rest mass density of the matter flelds (the 
Newtonian source term), but (among others) minus the square of the gradient of (post- 
Newtonian corrections) as well. Indeed, while the former is independent of c~^, the latter 
is proportional to and can naturally be interpreted as the energy-density (and/or 
the trace of the spatial stress tensor) of the gravitational fleld itself. Thus, to avoid 
confusion, we use the phrase 'post-Newtonian' in the following sense: (1) the spacetime 
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is static and asymptotically fiat, in which (2) when quantities are expanded as a series 
of then the zeroth order term are called Newtonian and the c~^^ order ones the A;th 
order post-Newtonian corrections. 

In the relativistically corrected Newtonian theory of gravity the volume integral of the 
source terms for the Newtonian potential gives a well defined (i.e. free of the ambiguities 
coming from the Galileo-Eotvos experiment) expression for the energy of the source plus 
gravity system even at a quasi-local level. Motivated by this observation, we show that 
in static spacetimes in exact general relativity the extra structures above are enough to 
introduce a notion of 'effective' quasi-local energy expression simply as the integral of 
the sum of all the source terms for 0. This energy is shown to tend at spatial infinity to 
the ADM energy as a monotonically decreasing set function. We give an explicit form of 
this 'effective' quasi-local energy in static, spherically symmetric spacetimes, and under 
certain energy conditions we prove its positivity and its rigidity, i.e. that its vanishing 
implies flatness. 

Since the idea behind the 'effective' quasi-local energy expression is exactly analogous 
to the quasi-local energy in the relativistically corrected Newtonian theory, we believe 
that the qualitative behaviour of the 'effective' quasi-local energy expression reflects some 
universality: we expect that any physically acceptable quasi-local energy expression in a 
static, asymptotically fiat spacetime should tend to the ADM energy as a decreasing set 
function at spatial infinity. 

In Section 2 we discuss the issue of energy both in Newtonian theory and in the 
relativistically corrected Newtonian theory. Then, in subsection 3.1, the Newtonian limit 
of Einstein theory is reviewed. We find that as a relativistic correction not only the 
energy density, but the spatial stress of the gravitational field also contributes to the 
effective source in the exact field equations for the Newtonian potential. We emphasize 
that this result is exact, i.e. no approximation is used. As far as we know this is a new 
post-Newtonian correction from GR, which has not been considered so far. Then the 
'effective' quasi-local energy is introduced and analyzed in subsections 3.2 and 3.3, and 
compared with the Misner-Sharp and Brown- York expressions in subsection 3.4. Section 
4 is devoted to the discussion of the potential implications for more general quasi-local 
energy-momentum expressions. 

The sign conventions for the metric and the curvature of [2U] are used. In particular, 
the signature of the spacetime metric is (+, — , — , — ), the curvature and Ricci tensors and 
the curvature scalar are defined by —R^'bcdX^v'^w^ := v^V c{w^V dX"-) — w'^Vc(f'^VdX") — 
[v,wY'VcX°', Rjyii '■= R°'bad and R := Ratg""^, respectively. Thus Einstein's equations take 
the form Gab '■= Rah ~ \Rgab = ~i^Tab — ^gabi whcrc A is the cosmological constant and 
K := SttG/c^ with Newton's gravitational constant G and the speed of light c; i.e. we use 
the traditional units. 

2 Gravitational energy in Newton's theory 
2.1 Newton's theory 

In a given inertial frame of reference the gravitational field is described by a scalar function 
(f) of the flat 3-space M^, for which the field equation is the Poisson equation 

-DaD^cj) = AnGp. 
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Here p : ^ [0, oo) is the rest mass density of the matter (source), and Da is the flat 
covariant derivative operator in the 3-space. Note that we use the negative definite fiat 
metric hab here, consistent with our conventions. If D C is any open subset with 
compact closure and a smooth boundary S = dD, then by the Gauss theorem and the 
field equation 

where denotes the outward pointing unit normal of S and dS is the induced area ele- 
ment. Thus the rest mass of the source in Newtonian theory of gravity can be rewritten 
into a 2-surface integral, like the charge in electrostatics. Following the analogy with elec- 
trostatics, we can introduce the energy density and the spatial stress of the gravitational 
field itself, respectively, by 

U--=^k-(DMD.^) = -^\DM\ (2.2) 

Sal, := ( (DA) (DA) - (DA) (D'A) )■ (2.3) 

In fact, the integral of f/ on S is just the work that we should do to form e.g. a spherical 
body by bringing particles together from infinity. Moreover, since gravitation is always at- 
tractive, in Newton's theory this is always a binding energy, and hence its sign is negative. 
The divergence of the stress tensor, together with the Poisson equation, yields the 'force 
density': -DqS"^ = ■^^{DaD^'cfyD'^cf) = —pD^cp. Note also that the average 'gravitational 
pressure' is just one-third of the gravitational energy density: 3P := —h'^^J^ab = U. 

However, by the Galileo-Eotvos experiment there is an important difference between 
electrostatics and gravitation. Namely, by this experiment the inertial and gravitational 
masses of the particles are strictly proportional to each other. Hence, in particular, the 
gravitational and inertial masses of the test particles drop out from the equations of 
motion in the gravitational field, yielding an ambiguity even in the notion of the gravi- 
tational force -De0: it is not possible, even in principle, to make a distinction between a 
uniform gravitational field and a uniform acceleration of the frame of reference. There- 
fore, at any given point of the 3-space the gravitational force De4> can be transformed to 
any given value, e.g. to zero, by an appropriate change of the reference frame. Thus the 
ambiguity in the gravitational force is De(f) ^ De(f) + ae, where is an arbitrary constant 
covector field in 3-space. It is only the second derivative DaDf,(j), the tidal force, that has 
direct physical meaning. Consequently, the gravitational energy density (12. 2p can also 
be transformed to any given non-positive value, e.g. to zero at any given point by an 
appropriate change of the frame of reference. Similarly, the spatial stress (12. 3 p is also 
vanishing at that given point. On the other hand, the gravitational energy density (as 
well as the spatial stress) can be transformed to zero on an extended, open subset of the 
3-space only if the gravitational field is uniform there. If, however, we have some extra 
information about the structure of the gravitational field, e.g. that it is the gravitational 
field of a localized source, then the ambiguity can be removed from the gravitational force 
Da4> by requiring its vanishing at infinity. 

To summarize, we see that even in Newtonian theory of gravity the gravitational 
energy and spatial stress cannot be localized to a point, just as a consequence of the 
Galileo-Eotvos experiment, and the rest mass of the source can also be written as a 
2-surface integral. 
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2.2 Two relativistic corrections to the source 

According to the special theory of relativity mass and energy are not independent con- 
cepts, and we should associate a mass distribution to any distribution of energy in 3-space. 
In particular, in addition to the mass distribution u/c^ of the internal energy density u of 
the matter field, we should associate a mass distribution U /c^ to the energy density f|2.2p 
of the gravitational field, too. However, according to the principle of equivalence, a con- 
sequence of the Galileo-Eotvos experiment, any mass distribution is a source of gravity, 
independently of the nature of the mass. Thus, in particular, both the internal energy 
density of the matter and the gravitational energy density contribute to the source of 
gravity. Therefore, the source term on the right hand side of the Poisson equation should 
be corrected, and the field equation for the relativistically corrected Newtonian theory of 
gravity could naturally be expected to be 



Note that the gravitational energy density reduces the magnitude of the source, because 
that is a binding type energy. As a consequence of (12. 4p we have that 



i.e. now it is the total energy of the source plus gravity system in a given domain D that can 
be rewritten into the form of a 2-surface integral. Note that while the gravitational energy 
density is ambiguous, this quasi-local expression for the energy of the matter plus gravity 
system is free of this ambiguity, just because the flux integral on S of any constant covector 
field ag is zero. This in itself already justifies the introduction and use of the quasi-local 
concept of energy in the study of gravitating systems. Moreover, in the source-free region 
(i.e. where p and u are vanishing) En is a decreasing set function, because then the 
integrand in the middle term is negative definite there. In particular, for a 2-sphere with 
radius r surrounding a localized spherically symmetric homogeneous source, the quasi- 
local energy is i?/),, = ^m(H-i^-|-^)-|-0(r~^), where, in terms of the rest mass M and the 
radius R of the source, the total mass parameteris m = ^(l-|^ + ^(^)^) + C'(c"®). 

3 Gravitational energy in static spacetimes 

3.1 The Newtonian hmit of Einstein's theory and two more rel- 
ativistic corrections 

In fT9], pp. 71-74, the Newtonian limit of Einstein's theory is defined through asymp- 
totically flat, static configurations. Thus let the spacetime be static, and K"^ be the (e.g. 
future pointing) timelike Killing field being orthogonal to the spacelike level sets t = const 
of a function t : M — t- M. These sets will be denoted by S^, and we write Ka = gVat 
for some function g on M. If := KaK"^, then t"" := f~^K"- is the future pointing unit 
timelike normal to the hypersurfaces S^, and hab '■= gab ~ tatb is the induced (negative 
definite) metric on S^. Let Da denote the corresponding intrinsic Levi-Civita covariant 
derivative operator on S^. Then by a straightforward calculation it is shown in |19] . pp 
72, that the length / of the Killing field satisfies the 'field equation' W^DaDbf = fRab^t^- 




(2.4) 




(2.5) 
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Next, let us define the energy density /i := Tabf^t^ of tlie matter fields seen by the 
static observers, decompose it into the sum of the rest mass and internal energy densities 
as /i = c^p + u, and write the trace of the energy- momentum tensor as Tabg"'^ = A* ~ 3p. 
(Thus —3p denotes the trace of the spatial stress tensor aab '■= PaPb'^cd with respect to 
the negative definite hab- Here, P^" denotes the obvious projection to S^.) Introducing the 
scalar field := (? In /, the field equation for / above, together with Einstein's equations, 
gives 

-/i'^''D,A0 = 47rG'p+^(n + 3p-^) +^/i'^''(Z}„0)(A0). (3.1) 

Comparing this equation with 02.41) we see that apparently we recovered 02. 4p with the 
relativistic correction term U /c^ dictated by the Galileo-Eotvos experiment (even with 
the correct sign), together with an additional relativistic correction: the trace of the 
spatial stress (as well as the cosmological constant) also contributes to the effective source 
(as already noted in [IH])- Nevertheless, the relative weight of the gravitational energy 
density term in the effective source is twice that of f l2.4p : the last term on the right hand 
side of 03. ip is rather than the expected ^^U. However, as we saw in subsection 

12. H we can associate with the Newtonian gravitational field not only energy density but 
also spatial stress, and the corresponding average pressure P is one-third of the energy 
density. Thus the 'extra' gravitational energy density in 03. ip can be written as 3P, i.e. 
the last term of 03. ip has the form ^^{U + 3P). Hence there is a fourth relativistic 
correction to (12. 4p : the gravitational stress also contributes to the effective source of 
gravity. Note that this correction is obtained in the exact theory, independently of any 
approximation method. (In addition, the intrinsic geometry {Htyhab), and hence the 
Laplacian —h°'''DaDi,, is not fiat. This can also be considered as an additional correction 
to 02. 4p . but it does not seem to be possible to formulate its deviation from the fiat-space 
Laplacian of in a gauge invariant way.) 

Another (and perhaps more direct) derivation of 03.11) could be based on the evo- 
lution parts P^P^{Gcd + i^Tcd + Xgcd) = of Einstein's equations in the standard 3+1 
decomposition. If C,"' = Nt"' + A^'* is an evolution vector field which is compatible with 
the foliation Sf, then, using the Hamiltonian constraint t°-t^{Gab + ^gab + i^Tat) = 0, the 
evolution equations are equivalent to 

Xab = N(^-^Rab + 2XacX^b - XXab) + L^Xab " DaD^N + 

+ XNhab + kN [-aab + ]^(y"chab + ^yU/iafe) , (3.2) 

where ^Rab is the Ricci tensor of the intrinsic 3-metric hab and L^^Xab is the Lie derivative 
of the extrinsic curvature Xab of along the shift vector field. Choosing the leaves 
of the foliation to be the hypersurfaces to which the Killing field K"^ is orthogonal the 
extrinsic curvature is vanishing, and choosing the evolution vector field to be the Killing 
field itself, the lapse will be the length of the Killing vector and the shift will be zero. 
Then 03. 2p takes the form 

- DaDbf = f [^Gab + H(Tab + ^^^(/U + ^P)hab^ , (3.3) 

where ^Gab denotes the Einstein tensor of the spatial metric hab and we used the Hamil- 
tonian constraint 
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= Kfx + \. (3.4) 

(Since in static spacetimes the local momentum density f^TabPc of the matter fields is 
vanishing, the momentum constraint is satisfied identically.) Taking the trace of fl3.3p 
and using (p instead of / we recover fl3.ip . In what follows we need the full fl3.3p rather 
than only its trace, and we consider (I3.3p -( l374|) to be the field equations rather than only 
(I3.ip . In fact, in the static case fl3.3p - (l3.4p are equivalent to Einstein's equations, the field 
equations for the Newtonian potential (p and the spatial metric hab- 

3.2 The 'effective' quasi-local energy for static configurations and 
its spatial infinity limit 

By the Galileo-Eotvos experiment any kind of energy is a source of gravitation. Thus, 
motivated by expression (12.51) of the relativistically corrected Newtonian theory, in exact 
general relativity in static spacetimes it is natural to define the total, 'effective' energy 
of the static matter+gravity system in a subset D G T^t, seen by the static observers 
as the integral of all the source terms on the right hand side of (13. ip : 

Note that contains not only the energy of the gravitational 'field' and (all kinds of) 
energy of the matter source, but the trace of the spatial stress of the source and the 
gravitational 'field' as well. Thus, apart from the cosmological term, the structure of the 
post-Newtonian part of the volume integral, u + U + 3{p + P), shows some resemblance 
to enthalpy rather than to the internal energy density. However, this combination seems 
to deviate from the standard form of enthalpy, too, which would have the structure 
u + U + {p + P). (For a more detailed discussion of the analogy of gravitational energy 
with the thermodynamical ones, see [S].) If the source is compactly supported in some 
Dq C C Et, and the cosmological constant is non- negative, then outside Dq the total 
gravitational energy is strictly decreasing with increasing domain D of integration. 

If the spacetime is asymptotically fiat (in which case A = 0) and the hypersurfaces 
extend to the spatial infinity, then Eo^, the quasi-local energy associated with the solid 
ball of radius R (or equivalently to the sphere Sr = ODr), tends to the ADM energy 
in the R ^ oo limit. To see this, it is enough to show that Ed^ tends to Komar's 
expression because it is known that the Komar expression built from the static Killing 
field K"^ (normalized to one at infinity) tends to the ADM energy (see e.g. [21] )• Recall 
that Komar's integral on a closed spacelike 2-surface S (with a timelike and spacelike 
unit normal, t"' and v"", respectively, and satisfying taV"" = 0) has the form 

ls[K'] := - (fv^^K'^l-Sabcd = - (f{v''t'ViaKb])dS, (3.6) 
1^ Js ^ 1^ JS 

where dS := \t'^v^eefcd is the induced area element on S. Using the Killing equation and 
the form Ka = exp(^)i:a of the Killing field, the integrand can be written as 

2vH^VyaKb] = 2v'^{VaKby - v^t\VaKb + VbKa) = 2 exp (^j t^^^^a^- (3.7) 
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Substituting this into (13.61) and taking into account that exp(^) — )■ 1 at infinity, we find 
that limn^oo Ed = hniij^^oo 1sr[K'']. 

3.3 Static spherically symmetric configurations: Explicit form 
and positivity 

Let the fine element of the spatial 3-metric be written as dh"^ = —e^^dr'^ — R^{d9'^ + 
sin^ OdcfP') for some functions R and a of r and regular at the origin r = with -R(O) = 0. 
Then the components of the outward pointing unit normal v"" of the Sr '■= {r = const} 

2- surfaces are = e~°'6l, and the corresponding extrinsic curvature is proportional to 
the induced 2-metric: Uab = -R^(5^(5^ + sin^ Here the prime denotes the 
derivative with respect to r. The curvature scalar of the spatial 3-metric is 

37^ = (l + 2i?i?'e-2°a' - 2RR!'e"^^ - (i?')'e-'") ; (3.8) 

while the curvature scalar of the intrinsic 2-metric on is ^TZ = 2/R^. (To avoid 
confusion, in this subsection the scalar curvatures are denoted by TZ.) 

To give an explicit form of the 'effective energy' in terms of the geometrical quantities 
defined on the 2-surface Sr-, we use not only the trace, but also the v°'v^ component of 
the evolution equation (13. 3p . The former is 

= D^D'^f + f(^-K{^^ + 3p) - a) = -e-"(e-"/')' - 2^6-^"/' + f(^-K{^^ + 3p) - a), 
while the latter is 

= V%'DaDtf + f(^^GabV%' + KaabV^v' - ^/^(/i + 3p)) = 

= e-"(e-"/')' + / - (f + ^TabV^v" - i«:(/. + 3p)) . 

Here, in the derivation of the second equation, we used the constraint (13.41) and the 
expression ^GabV°"v^ = |(^^ + ^ab^""^ ~ ^'^) for the f ^ component of the Einstein tensor 
of the 3-space in terms of the intrinsic and the extrinsic curvatures of Sr- Comparing 
these two and assuming that R' ^ 0, we obtain the expression of f~^f' in terms of R, a 
and TabV°'v^. However, this is essentially the integrand of the 2-surface integral in (13. 5p . 
v'^Df.cf) = c^e~"/~^/', yielding the explicit form of Er in spherically symmetric spacetimes: 

Er = ^^e'^[R'{KTabV%' - A) + 1 - (i^O'e"^"). (3.9) 

This formula can be simplified slightly if we choose r to be the areal coordinate, i.e. 
r = R, and hence R' = 1. (13. 9p gives Er algebraically in terms of the functions in the 

3- metric and the radial pressure Tabv"'v^ on the surface Sr of radius r. 

In the rest of this subsection we prove that this energy expression is non-negative if the 
matter fields satisfy the 'energy conditions' {hiTab + Xgab)t"'t^ > and {K,Tab + \gab)v"'v^ > 
0. Moreover, we show that, under the slightly stronger conditions {nTab + )^gab)t°'t^ > 
{uTab + \gab)v"'v^ > 0, the vanishing of Er implies uTab = —\gab and the fiatness of the 
Cauchy development of the ball with radius r. Clearly, {uTab + ^gab)v°'v^ > ensures the 
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non-negativity of the first term between the brackets in (13. 9p . To show that the second 
term dominates the third we use the constraint equation (13 ■4p . Substituting (13. 8p into 
fl3^ and multiplying by R'^R', we obtain {R{R'fe-'^°'y = R' -R^R'{k^i + X). Integrating 
this from zero to r and using {nTab + )^gab)t"'t^ > 0, we obtain 

i?(i?')^e-2" = i?- / {Kfi + \)R^R'dr < R; (3.10) 
Jo 

i.e. that {R'Y < e^°, and hence that > 0. (Here we used our previous assumption that 
R is strictly monotonically increasing, i.e. no minimal or maximal 2-surface is present.) 

Conversely, by the pointwise non-negativity of K,TabV°'v^ — A and of e° — (-R')^e~° from 
ii^r- = it follows their vanishing on the 2-surface Sr, i.e. in particular that R' = e". 
Substituting this to (13.101) we obtain that Kfi + X = on the whole 3-ball of radius r, and 
hence that R' = e° also on the whole 3-ball (and not only on its boundary). Substituting 
this into the line element it becomes fiat, i.e. the initial data set on the ball of radius r is 
the trivial one, and hence its Cauchy development in the spacetime is also flat. Finally, 
by the stronger energy condition, uTabf^t^ = —X implies nTabV^'v^ = A, i.e. nTab = —Xgab- 

3.4 Comparison with other round sphere expressions 

Since the Misner-Sharp energy appears as a mass expression in the study of equilibrium 
states of cold, spherically symmetric stars (see e.g. Appendix 1 of |19|). this became 
the more or less generally accepted definition of quasi-local energy on round spheres (i.e. 
on spherically symmetric 2-surfaces in spherically symmetric spacetimes). (Note that by 
spherical symmetry the spatial part of the energy-momentum is expected to be zero, and 
hence the mass is equal to energy.) In the line element of the previous subsection this takes 
the form Ems{^) = —r{l — e~^"), where r is the areal coordinate; while the Brown- York 
energy is EBYij") = ^'"(l ~ e~°). (N.B.: On round spheres the Bartnik, the Dougan- 
Mason and the Penrose masses and the Hawking, the Geroch and the Kijowski energies 
reduce to EMs{f)- On the other hand, the Brown- York energies (with all three choices 
for the reference configurations), the Epp, the Kijowski-Liu-Yau (which is Kijowski's free 
energy) and the Wang-Yau expressions reduce to EBvif)- For the details see e.g. [17].) 
Since the Hawking energy is a gauge invariant measure of the bending of the light rays 
orthogonal to the 2-surface (see subsection 6.1.1 of [IZ]), the Misner-Sharp energy can 
also be interpreted in this way. 

On the other hand, the 'effective' quasi-local energy was introduced as the integral of 
the effective source for the Newtonian potential seen by the static observers, i.e. there is 
a different concept of energy behind this: it is a measure of the effective source of gravity, 
including also the gravitational self-interaction. This yields that, in addition to the extra 
pressure and cosmological terms in (13. 9p . it has an extra overall weight function e" with 
respect to Ems{'^)- Since by the positivity proof above it is not less than 1, the energy 
Ej. is never less than the Misner-Sharp energy. 

In the Schwarzschild solution the Misner-Sharp energy is the constant — m for any 
r > 2ni, the 'effective' quasi- local energy is Er = ^—7=^, while the Brown- York energy 

V ^ r 

is EsYir) = ^^(1 — ^1 — ^). Thus, the Schwarzschild mass parameter m is only the 
'bare' mass, while in E^ this 'bare' mass is 'dressed' by the inverse local redshift factor 
and it tends to the ADM energy from above. Similarly, i?By(r) is also monotonically 
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decreasing. To see their more detailed asymptotic structure let us expand them as a 
power series of 1/r near infinity. We find 




(3.11) 



(3.12) 



which, by m ~ Gc^^ (see the last sentence of subsection 12. 2p . are the post-Newtonian 
expansions at the same time. Thus the Brown- York energy reproduces the Newtonian 
energy expression even in the order. On the other hand, by general relativity it 
is twice the Newtonian gravitational energy (more precisely the sum of the Newtonian 
gravitational energy and the average pressure) that appears as the source of gravity in 
the first post-Newtonian order (see (13. ip and the subsequent discussion). The factor 2 
in front of the post-Newtonian energy term in (I3.12p is simply a manifestation of the 
contribution of the gravitational spatial stress. Since, however, is a measure of the 
source seen by the static observers, it diverges at the horizon as it could be expected. Thus 
our 'effective' quasi-local energy may provide a physically reasonable notion of energy in 
the region where the Killing vector is timelike, i.e. outside the event horizon. 

4 Discussion 

In the literature there exist lists of a priori expectations on how a physically reasonable 
quasi-local mass or energy-momentum expression should behave |17[ [22], e.g. at spatial 
infinity, in the presence of spherical symmetry, or on the event horizon of black holes. 

One such additional natural expectation could be the compatibility with the results in 
the relativistically corrected Newtonian theory, where the quasi-locally defined energy of 
the matter plus gravity system tends to the total energy, measured at infinity, as a strictly 
decreasing set function. Thus, to be able to make this comparison, we need to define the 
Newtonian limit of general relativity. One of the several possible definitions is based on 
static spacetimes that are asymptotically fiat at spatial infinity |19| . Then we can expand 
every quantity and formula as a series of c~^, and while the zeroth order terms give 
the Newtonian approximation, the coefficients of c~'^^ are the kth order post-Newtonian 
corrections. In fact, since near spatial infinity in an asymptotically fiat spacetime the 
matter fields and the dynamics of both the matter and the geometry die off rapidly, the 
quasi-local energy-momentum/mass expressions could be expected to behave like in static 
asymptotically flat spacetimes. Though this notion of post-Newtonian approximation is 
more restrictive than the usual one (see [23]), the advantage of this is that the presence 
of the static Killing fleld provides extra geometric structures that make the subsequent 
analysis technically much easier and unique. In particular, they make it possible to deflne 
energy at the quasi-local level and to be able to compare energies (and not only masses) 
associated with different 2-surfaces unambiguously. 

Since in static spacetimes the gravitational contribution to the total energy is nega- 
tive definite, the quasi-local expressions should tend to the ADM expression as strictly 
decreasing functions. Therefore, in particular, if E{r) is any such gravitational energy 
expression evaluated on a large sphere of radius r near spatial inflnity and E^. denotes 
the 'effective' energy, also at r, then the asymptotic form of the former could be expected 
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e.g. to be E{r) = Er + 0(r~^) or E{r) = EBvif) + 0(r~^), depending on the concept of 
energy that is behind the actual notion E{r) of quasi-local energy. 

However, this new requirement is in conflict with two of the previous ones in [T7[ [22]. 
First, this contradicts to the expectation that for round spheres the quasi-local energy 
should reduce to the Misner-Sharp energy since the latter is increasing (or constant). 
Thus in static spherically symmetric configurations Er could be an alternative to the 
Misner-Sharp expression. 

Second, if the quasi-local mass should really tend to the ADM mass as a strictly 
decreasing set function near spatial infinity, then the Schwarzschild example shows that 
the quasi-local mass at the event horizon cannot be expected to be the irreducible mass. 
In fact, since both the ADM and irreducible masses are and the quasi-local mass must 
be strictly decreasing, there would have to be a closed 2-surface between the horizon and 
the spatial infinity on which the quasi-local mass would take its maximal value. However, 
it does not seem why such a (geometrically, and hence, physically) distinguished 2-surface 
should exist. 

This work was partially supported by the Hungarian Scientific Research Fund (OTKA) 
grant K67790 and by the Royal Society of New Zealand with Marsden grant UOO-09-022. 
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